On the physical origins of the negative index of refraction 
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The physical origins of negative refractive index are derived from a dilute microscopic model, pro- 
ducing a result that is generalized to the dense condensed phase limit. In particular, scattering from 
a thin sheet of electric and magnetic dipoles driven above resonance is used to form a fundamental 
description for negative refraction. Of practical significance, loss and dispersion are implicit in the 
microscopic model. While naturally occurring negative index materials are unavailable, ferromag- 
netic and ferroelectric materials provide device design opportunities. 
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In 1968, Veselago suggested that a material with a 
negative index of refraction would enjoy certain pecu- 
liar properties. Foremost among these was negative re- 
fraction, where a ray incident on the interface of such a 
material would refract on the same side of the normal 
rather than away from it pj. The subject lay dormant 
until 1999, when Pendry proposed designs for magnetic 
metamaterials |2| that were subsequently implemented 
by Shelby et at, along with a dispersed electric meta- 
material, to demonstrate negative refraction in 2001 
Despite the fact that no naturally occurring negative in- 
dex material is available, there has been a surge in inter- 
est, particularly with regard to the prospect of creating 
a perfect lens |^. 

In this Letter, we explore the physical origins of the 
negative index of refraction. The formalism developed 
here differs from previous macroscopic descriptions of 
negative refraction in that it is developed from a micro- 
scopic model consisting of an array of electric dipoles and 
magnetic dipoles undergoing precession due to an exter- 
nal static magnetic field. By considering the extension 
of the microscopic model to a macroscopic electromag- 
netic description, a physical basis for why the propaga- 
tion speed of light, v p , should appear negative when the 
permittivity and permeability are simultaneously nega- 
tive is established. This extends the macroscopic argu- 
ments originally formulated by Veselago [lj, and, for ex- 
ample, applied under the context of a Drude model 
to include a physical, microscopic basis. Finally, we use 
the microscopic model to suggest some possibilities for 
the use of ferromagnetic and ferroelectric materials to 
enhance the capabilities of negative index metamateri- 
als. 

To illustrate how negative refraction might arise, we 
apply Fermat's principle of least time to the case 
of negative refractive index materials and interpret it in 
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terms of least action. Consider the path of least action 
for a stream of photons in vacuum (n\ = 1) incident on a 
homogeneous, isotropic material with an index of refrac- 
tion n 2 . The most probable path is determined by the 
path of stationary phase, which corresponds to an ex- 
ternum in the spatial derivative of the total travel time 
through all possible paths. This is well established for 
the case of positive refraction. From the diagram in Fig- 
ure ID^a), the optical path length in vacuum from source 
A to the interface point O is c\t\ = AO = yj 'a 2 + y 2 , 
and in the semi-infinite material from O to B it is c^ti = 
OB = yjb 2 + (I — y) 2 . To find the extremum in the time 
of travel from A to B, we form d[ti + t>i)/dy — 0, with 
a and b fixed at arbitrary values, which upon substitut- 
ing for the optical path lengths gives n\yjcQ\J a 2 + y 2 = 
n 2 (l — y)/cQy/b 2 + (I — y) 2 , with I — y > 0. Recognizing 
the trigonometric relations sin(6*i) = yj \J a 2 + y 2 and 
sin(#2) = (I — y)l \Jb 2 + (I — y) 2 , this can be rewritten in 
the familiar form of Snell's law: nisin(0i) = r^sin^). 
If we postulate that the photons refract to the other 
side of the normal, as depicted in Figure HJb), then 
I ~ V < 0, which implies that the angle of refraction be- 
comes — 62- The extremum for this optical path gives 
niy/coy/a 2 + y 2 = -n 2 \l - y\/c y/b 2 + {I - y) 2 , which 
can only be satisfied if n 2 < 0. Since nisin(0i) > 
and sin(— 6*2) = — sin(#2)) Snell's law is found to be 
valid for both positive and negative refraction, provided 
we allow for the possibility of a negative index of re- 
fraction. The curvature becomes d 2 (ti + t2)/dy 2 = 
ni/[c (a 2 + y 2 f' 2 ]+n 2 /[c a {b 2 + (l- y ) 2 f/ 2 ]. Forn 2 > 0, 
the curvature is positive, indicating that Fermat's result 
indeed gives the minimum time (and distance). Interest- 
ingly, for n 2 < 0, the curvature is negative, indicating 
a maxima. In terms of least time this result does not 
make much sense, since the path for negative refraction 
corresponds a maximum, but in terms of least action, the 
path is acceptable upon recognizing that negative index 
materials are causal in an energy sense but not in a tem- 
poral sense, an issue we address further in relation to a 
microscopic model. 
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FIG. 1: Illustration for a) positive and b) negative refraction 
for a light ray incident on a material with \n 2 \ > |m|. In a) 
< n\ < n 2 and in b) < rij < \n 2 \ and n 2 < 0. 



To provide a mechanism for the apparent augmenta- 
tion of the arrival time of the photons, we employ a mi- 
croscopic model that replaces the slab with a thin sheet 
of electric and magnetic dipoles, find the fields due to 
each source, and then apply superposition. We assume a 
Lorentz electric dipole resonance and a Bloch magnetic 
resonance. 

The field radiated in the a;— direction by a thin sheet 
of microscopic electric dipoles is 



ET = 



2c, 
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(3) 



with the relative electric susceptibility given by 
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The disparity between the frequency dependent speed 
of light in a material and the constant speed of photons in 
a vacuum is resolved by recognizing that photons imping- 
ing on a medium can drive resonances in that medium, 
which may radiate and contribute to the total scattered 
field. When there is a difference in phase between the 
source and radiated field, the wavefronts at the detector 
will appear to be advanced or retarded with respect to 
the source, and it is on this basis that the concept of 
negative velocity is explained. 

To relate our microscopic model to the index of re- 
fraction, which is a macroscopic phenomenon, we employ 
an approach developed by Feynman, thereby extending 
the physical origins of the index of refraction to include 
negative index materials @. Our system consists of a 
monochromatic plane wave source of radial frequency u> 
traveling at the speed of light in vacuum, cq, in the in- 
direction, incident on a slab of material of thickness d, 
and a detector on the other side of the slab, sufficiently 
far away. If the slab were removed, then the observed 
travel time to the detector, a distance x from the source, 
would be commensurate with the speed of light in vac- 
uum. When the material is present and an observed delay 
or advance in arrival time is assumed, such that the ob- 
served speed is co/n, consider that the electric field at 
the detector is 



E d z = £ e^ d( "- 
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(1) 



which is the source field of initial strength Eq (neglect- 
ing amplitude scatter in the dilute limit) multiplied by 
a phase term. Feynman showed that for dilute media in 
which the phase term is small, Q can be written in a 
more suggestive manner as 
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Written this way, it is clear in that the plane wave at 
the detector is the sum of the source term (the incident 
field without the material) and a material response term 
(the scattered field). This representation is thus consis- 
tent with the Born approximation for scattered photons. 



where Xe = Nq 2 /me uj 2 is the static susceptibility, N is 
the oscillator density, q the dipole charge, e the permit- 
tivity of free space, m the reduced mass of the charges in 
the normal mode that results in radiation, ojq the elec- 
tric resonance frequency, and T is a phenomenological 
damping constant. 

Also within our slab is a sheet of magnetic dipoles 
with a static external field Hq applied in the i-direction. 
Equivalently, we could assume a slab of ordered spins, 
where the static magnetic field is provided by dipolar 
coupling, as in the case of a ferromagnetic material. The 
equation for the field radiated by a magnetic dipole sheet, 
E s z m , can be shown to be of the same form as (0 with 
the replacement of \e with the relative magnetic sus- 
ceptibility for a Bloch resonance, which is given by 



Xu 



Xmo r> t 



i + (w - n ) 2 T 2 



((fi - w)T 3 + i), 



(5) 



where x Mo is the ratio of the magnetization and the static 
magnetic field Hq that induced it, flo is the magnetic 
resonance frequency, and T2 is the spin dephasing time. 

The field due to both the radiating electric and mag- 
netic dipoles, using © and its magnetic analog, is then 



iujd 
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(Xc + X^Eoe^-^. 



(6) 



Comparing the microscopic equation of © with our pre- 
vious macroscopic expression in PJl. the index of refrac- 
tion can be expressed in terms of the microscopic suscep- 
tibilities as 



n = 1+ ^(Xe +XaO- 



(7) 



We will employ some foresight and recognize that (JJJ 
is the same as the index of refraction derived from the 
macroscopic Maxwell's equations for dense media of arbi- 
trary thickness in the limit of small susceptibilities. Ac- 
cordingly, we assume the usual form for the index of re- 
fraction, n = y^ErHr, where the relative permittivity is 
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e r = 1 + Xe and the relative permeability is \i r = 1 + x^, 
so as to not limit our discussion to dilute materials with 
weak susceptibilities. What we gain by doing so is that 
a layer of dipoles can augment the source field for subse- 
quent layers of dipoles. 

To describe the effect of the microscopic susceptibil- 
ity on wave propagation, we introduce the permittivity, 
e r = \e r \ exp(i(f) s ), and permeability, /Lt r = |/i r | exp(i^), 
in polar form into the index of refraction (n = \n\e 1 ^" — 
(|e r ||/v|) 1 /V W ' £+<W / 2 ), and wave admittance (Y = 
\Y\e** = y (kr|/|Mr|) 1/2 e^-^)/ 2 ). Y = (e /Mo) 1/2 
is the free space admittance, the inverse of the wave 
impedance. The domain for the index of refraction and 
the wave impedance come from the range of the permit- 
tivity and the permeability. Accordingly, with the loss 
implied by the damping in 10} and (J5J, the permittivity 
and permeability are restricted to the upper-half of the 
complex plane, in order for a propagating wave to lose 
energy (under the assumed e~ lut convention). Thus, (j) Etfl 
is in the interval [0, tt], and the index of refraction shares 
the same interval, since < (<f> e + (j>^)/2 < tt. However, 
the domain of the wave impedance phase is [— 7r/2, tt/2], 
since — n/2 < (<^> M — (f> e )/2 < n/2. Note that loss becomes 
more significant in the neighborhood of the resonances. 
However, even perturbational loss poses limitations on 
the potential amplification of evanescent fields in a neg- 
ative index slab |10| . 

The effect of the index of refraction on the source 
field as it propagates through the slab is seen by in- 
troducing the polar form of the index of refraction into 
and replacing the thickness d with the indepen- 
dent coordinate system in Figure ^ From , ignor- 
ing the time dependence, and assuming d — x, we get 

E^{x) = £; oe (- w:c M/ c o)sm0„ e i(wz|n|)/co)cos</>„^ gy uge f 

either the right hand rule or Faraday's law, we find that 
the magnetic field associated with propagation in the re- 
direction is H*(x) = -^{e 1 ^ e{(x). Thus, the wave 
admittance introduces a relative phase shift in the elec- 
tric and magnetic field, which has consequences on power 
flow. 

In the case of a material with only electric dipole cou- 
pling, the value of <p e is restricted to the interval [0, 7r] 
and </> M = 0. This limits the range of 4> n and <j)y to the 
interval [0, 7r/2]. However, if we allow for the simultane- 
ous existence of magnetic and electric dipole coupling, we 
find that the interval for the phase of the index of refrac- 
tion becomes [0,7r] and the interval for the wave admit- 
tance becomes [— 7r/2, 7r/2]. The most interesting result 
of this additional phase contribution comes from consid- 
eration of a material that without the magnetic dipole 
resonance would exhibit a bandgap for the frequency un- 
der consideration, assuming negligible loss, resulting in 
an evanescent wave. If we introduce a small magnetic 
phase contribution, then we find a return to a propa- 
gating wave solution. In the lossless electric dipole only 
case, the field does not propagate and does not carry real 
power because the electric and magnetic fields are 7r/2 
radians out of phase. However, when both electric and 



magnetic dipole coupling is present, the phase contribu- 
tion from the magnetic dipole can 'undo' the deleterious 
effects of the electric dipole by partially restoring the rel- 
ative phase of the electric and magnetic fields such that 
power can flow. In the remarkable case that <p £ and </> M 
are both equal to tt (or 0), it follows that <p n = tt (or 
0) and (f>Y — 0, and we find that propagation results 
with no attenuation, giving a time average Poynting vec- 
tor S = x±\E \ 2 Y ^/\e r \/\u r \. If the magnitudes of the 
relative permittivity and permeability are chosen to be 
identical, then unimpeded propagation takes place, i.e., 
there is no back-scattered wave. 

The negative refractive index range, <p n in [it/2, tt], cor- 
responds to the frequency range where (Xe+Xu)/2 < — 1, 
with Xq = x' q + i-Xq ■ These overlapping frequency bands 
for the resonances in |@J and © define the frequencies 
where propagation can occur. Within this band, a neg- 
ative phase velocity occurs, i.e., v p — co/n' < 0, where 
n = n' + in" . This means that with advancing time, the 
wave crests move in the — x-direction. However, in this 
same frequency range, the group velocity (v g ), which de- 
scribes the power flow of a wave packet, is positive. This 
can be established from v g = c/(n' + wdn'/dw) with use 
of J7J), (QJ, and ©. As a consequence, (uj/2)(dx' e /doj + 
dX^/dw) > |l+(Xe+xP/2| in the frequency range of neg- 
ative refractive index. Therefore, conservation of power 
holds, i.e., with incident power in the ^-direction, there is 
power flow in the x-direction in the negative refractive in- 
dex medium. Furthermore, causality can be established 
based on power flow, equivalent to the procedure used 
in dev elop ing the Kramers-Kronig relations for permit- 
tivity [111 . Note that conservation of energy is thus the 
basic metric for causality, and others such as least time 
based on phase (velocity), for example, Fermat's least 
time argument, require a revised interpretation from the 
standard view for positive index materials. 

The formalism developed here has direct application to 
ferroelectric and ferromagnetic resonances as either the 
electric dipole or the magnetic dipole source in a neg- 
ative index material. A good metric for the resonance 
frequency of ferromagnetics is the Bohr magneton di- 
vided by Planck's constant, which can conveniently be 
written as ~ 14 GHz/Tesla ! 9,f- This places the possi- 
ble range of resonance frequencies somewhere between 1 
and 100 GHz. Ferroelectrics do not have such a conve- 
nient metric, but the resonant frequencies tend to be in 
10's of THz Antiferromagnets can sometimes have 
resonance frequencies near a THz, but they have rel- 
ative susceptibilities similar to paramagnetic materials, 
which are only slightly greater than unity. It may, how- 
ever, still be possible to observe a negative index with an 
antiferromagnetic-ferroelectric composite. Within their 
respective frequency range, ferroelectrics and ferromag- 
nets offer very strong resonances, with susceptibilities ap- 
proaching up to 10, 000 in both. The frequency range 
over which the material may have a negative index is 
given by the extent of the bandgap in the limiting ma- 
terial, i.e., the one with the smallest bandgap, where the 
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bandgap is the region between the resonance frequency 
(wq) and it's conjugate frequency (a/) . In ferroelectrics, 
the latter is given by the Lyddane-Sachs- Teller relation, 
u>' 2 /u>q — So/eoo, where ojq is the resonance frequency 
usually associated with a transverse optical phonon in a 
ferroelectric crystal, eo is the static dielectric constant 
and £oo is the high frequency permittivity usually asso- 
ciated with electronic resonances |8(. In ferromagnetics, 
the bandgap is equal to f2o(Moo — 1), where /Zoo is the 
relative permeability above resonance 0]. These mate- 
rials offer an alternative to providing both electric and 
magnetic resonances with a metamaterial. Furthermore, 
ferromagnetic materials can be employed now with exist- 
ing metamaterials that operate normally around 10 GHz. 
One advantage in doing so is to exploit the modest tun- 
ability of the ferromagnetic resonance and the associated 
bandgap to provide a tunable negative index that could 
possibly be used in a device, e.g., a switch that refracts 
light one way or the other, based on the static magnetic 
field applied. Ferroelectrics also offer benefits for applica- 
tion to existing THz frequency split ring resonators 
in that the electric field of an electromagnetic wave can 
be imaged directly as it propagates through the ferroelec- 
tric material [14j , 

To summarize, by determining the path of stationary 
phase for a stream of photons incident on an arbitrary 
dispersive medium, we found that negative refraction oc- 
curs when the index of refraction is negative, indicating 
that the photons traverse this path with what seems to 
be a negative phase velocity. While Fermat's least time 
principle correctly predicts Snell's law for refraction at 
an interface with a negative refractive index material, it 
does not produce a least time solution, a consequence of 



the negative phase velocity. The microscopic model we 
have used for electric and magnetic dipole interactions 
provides a basis for the collective oscillator phase shift 
that results in negative phase velocity in the case of neg- 
ative refractive index. We found that when the electric 
and magnetic dipoles in this model have nearly the same 
resonant frequency, the augmentation of the electric com- 
ponent by the Lorentz dipoles is partially compensated 
for by the Bloch dipoles acting on the magnetic compo- 
nent, where the phase shifts in the electric and magnetic 
components of the field are directly related to the phase 
of the microscopic electric and magnetic dipoles. While 
any simultaneous resonance in the electric and magnetic 
constitutive parameters can provide a model for nega- 
tive refraction, the dilute microscopic interaction model 
presented, and its generalization to the dense limit, pro- 
vide a foundation for macroscopic interpretations. As 
a consequence, it is clear that dispersion and loss can- 
not be circumvented. Also, the dipole resonance models 
presented have a finite bandgap that is representative of 
physical systems such as ferroelectrics which can, with a 
concomitant magnetic resonance bandgap, provide nega- 
tive refractive index over a finite frequency range. While 
there are ferroelectric and ferromagnetic materials that 
may be candidates for negative refractive index applica- 
tions, achieving simultaneous bandgaps in an appropriate 
frequency range and having other satisfactory properties, 
such as low loss in the composite system, appears chal- 
lenging. On the other hand, use of one or the other phys- 
ical resonance, i.e., ferroelectric or ferromagnetic materi- 
als, in combination with a metamaterial implementation 
for the other resonance, is tractable and may yield prac- 
tical functionalities. 
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